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Abstract 

The light-like cusp anomalous dimension is a universal function in the anal¬ 
ysis of infrared divergences. In maximally {M = 4) supersymmetric Yang- 
Mills theory (SYM) in the planar limit, it is known, in principle, to all loop 
orders. The non-planar corrections are not known in any theory, with the 
first appearing at the four-loop order. The simplest quantity which contains 
this correction is the four-loop two-point form factor of the stress tensor mul- 
tiplet. This form factor was largely obtained in integrand form in a previous 
work for AA = 4 SYM, up to a free parameter. In this work, a reduction 
of the appearing integrals obtained by solving integration-by-parts (IBP) 
identities using a modified version of Reduze is reported. The form factor 
is shown to be independent of the remaining parameter at integrand level 
due to an intricate pattern of cancellations after IBP reduction. Moreover, 
two of the integral topologies vanish after reduction. The appearing master 
integrals are cross-checked using independent algebraic-geometry techniques 
explored in the Mint package. The latter results provide the basis of mas¬ 
ter integrals applicable to generic form factors, including those in Quantum 
Chromodynamics. Discrepancies between explicitly solving the IBP rela¬ 
tions and the MINT approach are highlighted. Remaining bottlenecks to 
completing the computation of the four-loop non-planar cusp anomalous 
dimension in AA = 4 SYM and beyond are identified. 
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1. Introduction 

Maximally (AA = 4) supersymmetric Yang-Mills theory (SYM) offers 
perhaps the best chance of truly solving an interacting four-dimensional 
quantum held theory. In addition, it is a proven stepping stone to pio¬ 
neer new computational technology, as the large amount of supersymmetry 
renders many computations much simpler than their non-super symmetric 
counterparts. In many cases, the resulting technology has transformed com¬ 
putational power in generic quantum held theories. For instance, in QCD, 
recent years have seen a boom of next-to-leading-order (NLO) computa¬ 
tions considered unfeasible only a few years before, inspired ultimately by 
Witten’s twistor string [l|. 

The focus of this paper is the computation of the Sudakov form factor, 
which is an observable that involves one gauge-invariant operator from the 
stress tensor multiplet and two on-shell massless states, 

T2 = {PI,P2\0\0). (1) 

In A/" = 4 SYM, this form factor was hrst discussed and computed at the 
two-loop order in 0]. It is noteworthy that the three-loop form factor was 
hrst computed in QCD j^, almost three years after the master integrals 
were found in j^. This result was then adapted to provide the three-loop 
answer in A/" = 4 SYM in j^. The latter computation cleanly shows that 
the computation within A/" = 4 SYM is technically much more straightfor¬ 
ward, especially when employing modern unitarity-based methods instead 
of Feynman graph techniques. 

Since there is only a single scale in the problem of computing the Sudakov 
form factor, the dependence on this scale is hxed by dimensional analysis. 
Hence, the form factor evaluates to a function which only depends on the 
dimensional regularization parameter e (D = 4 — 2e), the coupling constant 
and the number of colors Nc- Moreover, the dependence on e is governed 
to a large extent by the known structure of infrared (IR) divergences. The 
divergent structure through for instance, is determined by the so-called 
cusp (or soft) anomalous dimension, which gets its name from its appearance 
in the computation of the light-like cusped Wilson line 0,0. The cusp 
anomalous dimension is a function which is universal for a given quantum 
field theory. 
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In AA = 4 SYM, an integral equation which determines the leading pla¬ 
nar part of the cusp anomalous dimension in principle exactly was derived 
from integrability [81]. This is to date the most powerful and precise mani¬ 
festation of the cluster of ideas known as the AdS/CFT correspondence . 
However, much less is known about the non-planar part of the cusp anoma¬ 
lous dimension. At weak coupling, an immediate problem is that the first 
non-planar correction appears at four loops and has never been calculated 
to date, in any theory. Further motivation comes from a conjecture 
that this correction may in fact vanish. This is based on the assumed com¬ 
pleteness of the dipole contributions, which is consistent with the analysis 
of IR divergences in Yang-Mills theories explored through three-loop order, 
see also 11[. However, there is evidence from the Regge limit analysis in 
0 that this naive dipole summation is not sufficient at four loops. It is 
important to settle this issue by explicit computation. 

The integrand of the form factor was obtained in 0 using color-kinematic 
duality IJ, [l5j, [16| as an ansatz generator. The coefficients in the ansatz 
were fixed by a choice of unitarity cuts, up to one remaining parameter. To 
evaluate the integral, a major next step is to perform the integral reduction, 
in particular using integration-by-parts (IBP) reduction 17], [l8|. This turns 
out to be highly non-trivial due to the complexity of four-loop integrals. The 
success of this reduction is a major achievement of this paper. In particular, 
this reduces the challenging non-planar cusp anomalous dimension problem 
down to the computation of an explicit set of master integrals. 

Based on the IBP reduction result, surprisingly, the free parameter in 
the integrand turns out to be truly free: it drops out of the full form factor 
after the reduction. This completes the determination of the integrand of 
the M = A SYM form factor at four loops. In other words, we prove that 
there is a one-parameter family of four-loop form factor integrands which 
satisfy color-kinematic duality. 

The master integrals are also explored independently by algebraic tech¬ 
niques. While the results are consistent with the IBP reduction in majority, 
some rare examples of mismatch present interesting open questions to under¬ 
stand further. The results from algebraic methods, with some qualifications, 
provide the master integrals of four-loop form factors in generic quantum 
field theory, including QCD. 

The remaining problem of computing the four-loop form factor in A/" = 4 
SYM is the evaluation of the master integrals. Many techniques are available 
to compute Feynman integrals, see 19| and references cited therein for an 


overview. However, the basis integrals of the four-loop form factor are rather 
non-trivial. They include non-planar four-loop integrals with up to quadratic 
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irreducible numerators, the majority of which, to our knowledge, have never 
been integrated. We present partial numerical results and briefly survey the 
challenges involved in completing the computation, while leaving the full 
evaluation to future work. 

This article is structured as follows. First, the connection between form 
factor and cusp anomalous dimension is reviewed in section [2j This sec¬ 
tion also introduces some of the language to study Feynman integrals and 
discusses the issue of graph symmetry. In section [3l we report on the IBP 
reduction of the Feynman integrals in the four-loop form factor in AA = 4 


SYM based on the Reduze package 20(]. The independence of the sum on 


the free parameter reported in 0 is elucidated. We then proceed, in section 
m to study the basis of master integrals using algebraic techniques explored 
in . A brief survey of the numerical evaluation is given in section [H Ap¬ 
pendices contain information on a three-loop cross-check and an alternative 
choice of basis integrals. In the supplementary material section, there are 
several ancillary files which contain results for the basis of master integrals. 


2. Review^ and setup 

In the following, extensive use will be made of the results in for the 
four-loop form factor in AA = 4 SYM, in particular of the graphics and tables 
in section 5 of that paper. The numbering of the integral topologies refers 
to that paper. 

In the course of writing this paper, a typo was discovered in Table 5 
of [l^: the entry for the color factor of integral topology 26 should be 
2 AN^ 6 aia 2 - III particular, this integral topology does not contribute to the 
planar cusp anomalous dimension. 


2.1. Relation between form factor and cusp anomalous dimension 

The Sudakov form factor is an observable which here involves one gauge- 
invariant operator from the stress tensor multiplet and two on-shell massless 
states, 

A-2 = (pi,P2|O|0). (2) 


In the present paper, we focus on the four-loop Sudakov form factor in AA = 4 
SYM, for which the integrand was obtained in 00 


^In Af = 4 SYM, form factors were first studied thirty years aRO in and revived 
in the past few years at weak coupling [1, [H, 0, [13, [H, |3, [13, HI, 113, 0, Hll, HI] and 
at strong coupling [l^ . [13,113] . See also recent developments concerning form factors of 
non-Bogomolnyi-Prasad-Sommerfield (BPS) operators [13. IItI. 11^. l39l| . 
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The Sudakov form factor for the stress tensor multiplet can be written 
as the tree-level form factor times a scalar function, 


J 2 F2(5,A^c). 


( 3 ) 


It is this function F 2 which will be computed. In momentum space, F 2 is 
a function of the only scale in the problem, = (pi +P 2 )‘^, the coupling 
constants, group theory Casimirs and e. Here, q is the momentum associated 
with the gauge-invariant operator. Note that it does not obey a mass-shell 
condition, in contrast to the momenta of the gluons, Pi=P 2 = 0 . 

The cusp anomalous dimension 7 cusp is related to the two-point form 
factor through the universal exponentiation of IR divergences 40|, [dlj, [42, l43| , 
which, in dimensional regularization, takes the form 


Log[T2] = 


^2\-le ■ 


-g^hcnsp 

4(/e)2 


+ 0{e-^) 


( 4 ) 


The cusp anomalous dimension can be calculated as a perturbative series in 
the coupling constant. Displaying the dependence on group theory quanti¬ 
ties up to fourth order, it reads 


Tcusp — 


E 


9 Tcusp 


= aig^CA+a2g'^C‘i+a3g^C\+g^ {cl\C\ + dii)+0{g^^) , 

( 5 ) 



00 with 


where Ca is the Casimir which, for the SU(iVc) gauge group, simply has value 
Nc- The quantity ^44 is a particular group theory Casimir invariant which, 
for SU(iVc), has the value In previous works, the constant s aj 

02 , 03 and were calculated in perturbation theory 
and using integrability [1]. In ’t Hooft’s planar limit 
g'^Nc fixed, it is clear that the first subleading color correction to the cusp 
anomalous dimension occurs at four loops. 

From the general exponentiation formula in equation (j4]), this non-planar 
correction can be isolated from the non-planar part of the form-factor at 
four-loop order at order ^ in the e expansion. Since the integrals present 
in the form factor naively diverge as this gives an important check on 


^The three-loop cusp anomalous dimension in A/” = 4 SYM was first obtained 
turally) in (48| by extracting the leading transcendentality part of the QCD result 
In this case, the cusp anomalous dimension is related to the anomalous dimension of the 
twist-two operator in the limit of spin j ^ 00 . 
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any computation of the non-planar form factor: the first six orders in the ^ 
expansion must vanish. 

The planar correction to the cusp anomalous dimension at four loops 
can be expressed as a sum of rational coefficients times either C(3)^ or tt® 
471, which are both constants of transcendentality six. Based on the 


maximal-transcendentality principle [481], it is natural to conjecture that the 
four-loop non-planar correction 04 ^ is given by a combination of the same 
transcendental numbers, if not vanishing. 


2.2. Feynman integrals 

An L-loop Feynman integral with n so-called indices ai,..., is an 
integral which can be written as 

/(ai, ...,an)^ I d^h ■ ■ ■ d^lL (l/Dir • • • , ( 6 ) 

where Di are inverse propagators. For future reference, we define the ab¬ 
solute value of the sum over positive and negative indices to be t and s 
respectively. The parameter s is the numerator power and t is the propaga¬ 
tor power. In this article, the indices are strictly integer valued. Integrals 
which have the same non-zero indices for propagators are said to belong to 
the same sector. In the special case where all propagator powers are one, 
the integral is referred to as the corner integral of the sector. 

The problem under study in this article consists of a series of four-loop 
twelve-propagator integrals. For each integral topology, one can construct a 
basis by starting with a parametrization of the loop momenta of the graphs. 
Then, one needs to pick six additional propagators to get a complete set. 
These additional propagators are important for expressing irreducible nu¬ 
merators. Given a numerator, one can express this uniquely into the basis. 
The choice of numerators is by nature somewhat arbitrary. If one focuses 
solely on additional propagators which are products of differences of two 
momenta, then this can be done by trial and error. A more systematic 
approach is needed, if the aim is to include graph symmetry. 

As an example, let us consider integral topology 26. Its propagators can 
be parametrized as: 

ll ll ll , ll , i-k + h? , {h - h)\ {-h+Pl)\ {-h+P2)\ 

{—k+h+k)'^ 1 (~^3+Pi+P2)^) {k—k+k~P 2 )‘^i {—k+k—k~kFpi+P 2 )‘^ ■ 

(7) 
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This topology has two independent permutation symmetries, as will be dis¬ 
cussed later. It can be checked that the following six additional propagators 
parametrize all irreducible numerators: 

{h — h)'^ , — ih — , {h—pif, ih—P2 )‘^, {h—P2f- (8) 


In section pVppendix Ct choices of propagators are explicitly given for all 34 
topologies. 

A useful step toward computing the Feynman integrals is the so-called 
a representation, see e.g. 19(]. The result is an integral over a parameters, 
one for each propagator. 


I{(ll , ■ ■ ■ , Un) CX 


1 


n.r[o. 


f 






-\-rj-Dl2-\FIU 


(9) 


where the normalization constant, which is unimportant for our purposes, 
has been suppressed. U and F are certain polynomials of the a parameters 
of homogeneity L and L-|-1, respectively. See [l9| for further discussion and 
for a Mathematica code to compute them from given propagators. As is 
customary, the labels on the Oj parameters correspond to the position in a 
given ordered list of propagators, such as the one in equation ([7]) for integral 
topology 26. Equation ([9]) can cover negative indices, that is, numerators as 
well. 


2.3. Graph symmetry 

Several of the graphs in the set under study have a graph symmetry: 
the graphs are mapped to themselves under a permutation of some of the 
edges and external legs. Exploiting graph symmetry can be of great help in 
simplifymg and cross-checking computations if it is present. 

In [131], graph symmetry was indicated briefly by the symmetry factor. 
These symmetries can be checked or found by using the built-in capabilities 
of Mathematica. These permutations can include permutations of the exter¬ 
nal legs. An alternative strategy is to compute the U and F polynomials of 
the a parameter integral, equation Then, one can simply check all pos¬ 
sible permutations and see if they leave these polynomials invariant. Due to 
the factorial growth of this problem, the feasibility of this analysis strongly 
depends on the algorithm: a naive implementation ran out of steam for ten 
propagators. A faster algorithm is presented in (53|. 

A symmetry of the graph of an integral can correspond to a symmetry 
of a corresponding Feynman integral. For a given parametrization of the 
loop momenta of the graph, a permutation symmetry of the graph may be 
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undone by a shift of the loop momenta, combined with a permutation of 
the external legs. If the latter leaves the integral invariant, the result is a 
symmetry of the Feynman integral. 

For the example of topology 26, the two independent permutations are 
generated from 


Iq 0/5, 

I4 —>■ —I4 + /s + /e ) 
h -^-k+Pi+P2, 
Pi ^P2, 


and 


le ~h + Pi ) 

I 5 + P2 , 

Ia I 3 ~ h 1 

Pi ^ P2, 


( 10 ) 


as can be checked explicitly in equation ([7|. Since the integrals, by dimen¬ 
sional analysis, only depend on (pi +P 2 )'^, they are left invariant by the 
permutation of pi and p 2 - 

The momentum map given above corresponds to a simple permutation 
pattern for a given list of propagators. Let us consider, for instance, the 
scalar integral in topology 26, with the propagators as listed in equation 
©• The two independent permutation symmetries of the graph correspond 
to the following cycles: 


{{1,7}, {2,8}, {3,6}, {5,11}, {9,12}}, 
{{1,2}, {3, 9}, {4,10}, {6,12}, {7,8}}. 


Together with the trivial permutation and their product, these form a four- 
element representation of the permutation group. A given propagator can 
be mapped to other propagators in the set, but only if they are in the orbit 
of the given permutation. The orbits of the above permutation group are 


{{1,2,7,8},{3,6,9,12}, {4,10}, {5,11}} . (12) 

This can be used to simplify the sector decomposition method employed in 
FIESTA, for instance, as will be shown in section EH 

Given a permutation symmetry of a particular topology, ideally, addi¬ 
tional propagators would be found which yield a complete basis of propa¬ 
gators and also incorporate the permutation symmetry. One can search for 
symmetric numerators systematically. Permutation symmetries either leave 
propagators invariant or interchange two propagators. Hence, the six prop¬ 
agators to be added to the set contain either zero, one, two or three pairs of 
interchanges and, consequently, six, four, two or zero invariant propagators. 
Starting from a generic linear combination of all four loop momenta and two 
external momenta, one can construct the most general pair combination as 
well as the most general invariant propagator in terms of a number of free 
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variables. Since propagators are only determined up to a sign by the appear¬ 
ance of the square, there are two times two possible (pairs of) polynomials 
which come out of this. Hence, all possible appearances of the permutation 
symmetry can be parametrized. 

Then, one checks if the resulting set of eighteen propagators is linearly 
independent. In several cases (e.g. topology 20), no set exists at all. In 
addition, there are integrals with multiple permutation symmetries which 
do not permit a choice of six additional propagators with all permutation 
symmetries manifest. In cases in which there are sets of numerators with 
explicit permutation symmetries, typically many free variables exist. These 
may then be fixed by hand aiming for propagators as simple as possible. 
Care should be exercised not to choose parameters such that accidental 
linear dependencies in the total set of propagators are introduced. 

For the example of topology 26, the momentum maps in equation (jlOl) 
do not map the additional numerators in equation ([8]) into a permutation of 
themselves. The following set of additional numerators has explicit permu¬ 
tation symmetry: 


(^3 -I- 2^5 — ‘ipif , 

{—h + 2^6 + P2)^ , 

(^3 — 2^6 — Pi + 2p2)^ , 


(/a -I- 2/5 — p2f , 

[h — h + h— P 2 f , (13) 

{h -k-pi +P2f ■ 


This example shows a general principle: the price of manifest permutation 
symmetry can be a much more complicated numerator structure. It would 
be very interesting to fix guiding principles for choosing sets of numerators. 

The considerations above result in two sets of propagators. For each in¬ 
tegral topology, there is a choice of eighteen propagators which can be used 
to express any integral in this particular topology. In particular, any numer¬ 
ator can be expressed in terms of the basis. In the first set, the aim was to 
include numerators as simple as possible, while, in the second, permutation 
symmetries were taken into account as much as possible. Obtaining the 
symmetric set is the result of a fairly lengthy computation. The simplest 
choice in the hrst set is listed in section (Appendix C). All results reported 
here were obtained with this set, unless explicitly mentioned otherwise. 


3. IBP reduction and its output 

The set of Feynman integrals in equation ([6]) is over-complete, as there 
are relations between different Feynman integrals. A particular example are 
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the IBP relations 0,0 , which follow from 


0 = / d^h...d^lL -^X. 


(14) 


Working out the right-hand side leads to a linear relation of different Feyn¬ 
man integrals. By solving a system of such equations, one may express a 
general Feynman integral in terms of some basic integrals. 

A different way of phrasing the problem is to envision the system of 
IBP relations as a large matrix equation, with the integrals combined into 
a vector. The standard way of solving a problem of this type is Gaussian 
elimination. However, the output of this algorithm for a non-invertible ma¬ 
trix problem, such as the one under study here, depends on the ordering 
of the integrals in the vector. This is the essence of Laporta’s algorithm 


for IBP reduction 5^: one picks an ordering of the integrals. This or¬ 


dering should be such that complicated integrals are expressed in terms of 
simpler ones in general. For instance, integrals with smaller values of the 
parameter s (measuring numerator powers) are preferred, as are integrals 
with less numbers of different propagators. In addition, smaller values of 
the parameter t (measuring propagator powers) are to be preferred. Within 
these general choices, in practice, the exact criteria can differ between two 
different implementations of Laporta’s algorithm. The results of these two 
implementations are, therefore, in general not the same, but related by a 
change of basis. 

Given a system of equations, the output of Laporta’s algorithm is a 
reduction of all integrals in a given set down to integrals which cannot be 
further reduced from the given equations using the given ordering. The left¬ 
over integrals are known as master integrals. These integrals depend on the 
given set of equations, although typically, if the given set of equations is 
large enough, the set of master integrals tends to converge. 

Various private and public imp lementations of Laporta’s algorithm exist, 
such as AIR M], FIRE M, 5^ and Reduze [^, 59]. See LiteRed 61] 
for an alternative approach to IBP reduction. We explored FIRE, Reduze and 
LiteRed in some detail for the four-loop form factor problem. Only Reduze 
succeeded in solving the problem, probably due to its better parallelisation 
implementation. 


3.1. Relevant implementation details 

Reduze is designed to run in massive parallel mode on a cluster. As an 
input, it takes a family of integrals. This family is treated as an ordered 


10 












set: lower integrals in the set are considered simpler. Integrals are reduced 
as far as possible to lower integrals. This cuts down on the computational 
workload considerably, especially in a large-scale problem such as the one 
under study here. As parameters, Reduze takes a range of propagator powers 
(the t variable) as well as a range of numerator powers (the s variable) and 
constructs out of these all possible IBP relations for a given sector. These 
are then solved sector by sector, starting from the simplest sector, which is 
the bottom sector of the first member of the integral family. The number of 
relations tends to grow very quickly along either the t or the s axis. 

A problem we encountered in the public version of Reduze in parallel 
mode was that it tended to crash when handling many large files on the 
cluster file system during the identity generation stage. Due to this prob¬ 
lem it is unfeasible to obtain an IBP reduction with the public version of 
Reduze, and a fix is needed. The problem was that a file was reported as 
not present in the file system, while a manual check thereafter did uncover 
this particular file. This is due to the internal file handling structure of 
Reduze H where all processes are allowed to read from and write to disc. As 
a work-around, the program can be forcecfl to sleep for a time-out of fifty 
seconds, followed by a 500-second time-out if the file is still not found. This 
two-stage timeout resolved the crashing problem due to disk space usage. 
The loss of productivity for a few minutes is a small price to pay to prevent 
the total collapse of the computation. 

During the running of the four-loop reduction, Reduze still crashed oc¬ 
casionally when handling integrals in sectors with many, typically twelve 
or eleven propagators. These integrals require a large memory usage each, 
which can easily overwhelm one of the nodes of the cluster. Sometimes this 
can be due to bad scheduling of processes, with all root processes which con¬ 
sume most memory running on a single node; this situation can somewhat 
be avoided at start-up by passing scheduling instructions to the message¬ 
passing-interface (MPI) program^ By contrast, for lower propagator count, 
the limiting factor for speed is mostly CPU time. These two situations re¬ 
quire opposite numbers of allowed parallel processes, which in turn requires 
occasional input by the user during running. It should be noted that, al¬ 
though Reduze works in parallel, there is a limit to the number of cores 
which are assigned to a single process while still increasing computational 


®We thank A. von Manteuffel for explaining this. 

^This can be done in the function set_job_status in the file jobqueue. cpp. 
®One uses -bynode in openmpi. 
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speed. This saturation is sector dependent, but certainly under a hundred 
cores. 

The chosen IBP relations are in the range t < 12 and 0 < s < 2, with 
some extensions to t = 13 for specific integral topologies. Since Reduze 
features a choice of used relations, there can be problems with the reduction 
at the outer edges of these relations: there can be unresolved integrals. A 
simple check of this is to observe the file size of the reductions in a sector 
and compare it to the average for other integral topologies in a sector with 
the same numbers of propagators. A sudden increase in file size tends to 
indicate unreduced integrals: the unreduced integrals appear with much 
more complicated prefactors than the others, typically involving polynomial 
ratios with polynomial orders an order of magnitude above the norm. This 
is more than just a nuisance. The abnormal prefactors tend to be highly 
divergent in the limit e —>■ 0, typically oc ^ or worse. Since the form factor 
itself is expected to diverge at worst as this implies massive cancellations 
among the integrals. This particular problem occurred in the basis of section 


Appendix C in integral topologies 20 and 26, for integrals with at most ten 
propagators. This problem was resolved using consistency conditions, see 
two paragraphs down. 


We reduced using the set of numerators in section Appendix C A sec¬ 


ondary computation using numerators with explicit permutation symmetries 
failed; at a certain eleven-propagator integral, the MPI protocol ran into an 
issue with maximal size of a single message, which is set by the maximum 
size of an integer on a 32-bit machine, see 6^ for an explanation of the 


problem and a possible resolution. This is a result of a choice of protocol 
and not related to the physical architecture used. In principle, this problem 
could be circumvented by rewriting the Reduze code to handle messages of 
this size. 

Since the set of numerators with the simplest structures was used, inte¬ 
gral relations from permutation symmetries are not explicit. This allowed 
a simple check on the results: take a graph symmetry relation and use the 
IBP reduction. If Reduze found the relation, the result must be zero. We 
checked permutation symmetry on maximal-propagator integrals with up- 
to-quadratic numerators. Surprisingly, one single additional relation was 
found which involved exactly all master integrals with abnormal prefactors. 
Plugging this relation into the reduction eliminated the abnormal prefactors, 
reducing them to much more normal-looking ones. 
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3.2. Results 


The computation led to a result after several months of computation on 
a fairly large cluster]^ The complete list of the master integrals obtained 
with this reduction is attached, both in the planar and non-planar sectors, as 
well as a sample result for topology 25. Overall, the integrals which appear 
after reduction are simpler. There are two integrals with two numerator 
powers. One of these contains only eight propagators and is easy to integrate 
numerically. The other, however, has the full twelve propagators and occurs 
in integral topology 26. It turns out that this particular integral can be 
integrated numerically (see section EH). 

Any basis of master integrals listed is arbitrary to some extent. For a 
given sector, an obvious first choice for a master integral is its corner integral. 
Beyond this, there is a choice of a (typically single) numerator or doubled- 
up propagator. The statistics in Table [U relates to the direct output of 
Reduze, taking into account the one extra relation found from permutation 
symmetry. The statistics is split into those integrals which contribute to the 
planar and non-planar parts of the form factor. 

Table 1: Master integral statistics of obtained IBP reduction. 


(a) planar form factor (b) non-planar form factor 


# props 

s = 0 

S=1 

s = 2 

# props 

s = 0 

S = 1 

s = 2 

12 

8 

6 

0 

12 

10 

10 

1 

11 

18 

2 

0 

11 

13 

3 

0 

10 

43 

9 

0 

10 

34 

10 

0 

9 

49 

1 

0 

9 

29 

1 

0 

8 

51 

4 

1 

8 

32 

3 

1 

7 

25 

0 

0 

7 

13 

0 

0 

6 

8 

0 

0 

6 

7 

0 

0 

5 

0 

0 

0 

5 

1 

0 

0 

sum 

203 

22 

1 

sum 

139 

27 

2 


In these tables, there is quite some double counting of integrals which 
appear both in the planar and non-planar sectors. Using the known result 
for the sum of the integrals in the planar sector, one can eliminate several of 
the integrals in the non-planar sector. This is more efficient than might be 


^Memory usage on single sectors can exceed 200 GB. More detailed information is 
available on request. 
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expected, as there are four integral topologies (21, 25, 30 and 31) which each 
appear with the same overall color factor, 2N^ + 24A(?. This can be used 
to eliminate these integral topologies from the non-planar sector in favor of 
the integrals in the purely planar sector and the known planar-form-factor 
result. The resulting sum contains more integrals (260 in total) than the 
non-planar sector by itself. 

The choice of master integrals in the provided results is taken from the 
output of Reduze. It is possible to choose other sets of master integrals, in 
particular the ones with numerators, to aim at simplifying, for instance, the 
prefactors of the integrals in the form factor or to make use of a permutation 
symmetry of the scalar skeleton of the integral. In particular, one aim could 
be to make these as weakly divergent as possible in the limit e —)■ 0. The 
power of divergence in this limit determines the order to which the integral 
is to be expanded to get the cusp anomalous dimension. 

Cross-check: Multiple reductions 

The IBP reduction was run twice, with a second reduction involving 
a reshuffling of the integral topologies to put planar topologies first. It 
was checked that the two reductions commute: reducing with the second 
reduction followed by the first gives the same result as reducing with the 
first, after taking into account the extra relation found from permutation 
symmetries. Before taking into account this relation, there are seeming 
mismatches between several of the integrals which appear in the full four- 
loop form factor: these mismatches are all proportional to the extra relation. 
This is a consistency check for the IBP reduction. 

3.3. Cancellations in the form factor after IBP reduction 

The IBP reduction displays two separate patterns of interesting cancella¬ 
tions for the form factor result: one pattern takes place within the separate 
integral topologies, while the second is between different integral topologies. 

The first pattern is that, for integral topologies 18 and 20, the contribu¬ 
tions to the form factor vanish after IBP reduction. Generically, the Feyn¬ 
man integrals in these topologies do not vanish; there are twelve-propagator 
master integrals. However, the specific combination of the integrals in the 
TV" = 4 SYM form factor reduce to zero. Considering the complicated nu¬ 
merators for topologies 18 and 20, this is rather nontrivial and not obvious 
at all at the integrand level. It should be noted that these integrals already 
do not contribute to the form factor: their color factors are also zero. Such 
cancellations, however, do not happen for the other four zero-color-factor 
topologies, 8, 11, 15 and 16, whose integrals reduce to nontrivial results. 
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It would be very interesting to find an explanation for these vanishing in¬ 
tegrals, as, in physics, a zero almost always has a physical and important 
symmetry explanation. 

Another intriguing pattern emerges when tracking the dependence on 
the single free parameter in the form factor integrand which was left in 
[l3 |. There, the free parameter in the form factor multiplies a combination 
of 24 out of 34 different integral topologies which obeys all the physical 
constraints imposed in separately. In particular, this combination is 
color-kinematic dual by itself. One may expect that, by applying other, 
more complicated cuts which are not exploited in [l^ . this parameter may 
be fixed uniquely. The surprise is that, when summed with the correct color 
and symmetry factors (and only then), the dependence on the free parameter 
drops out of the form factor after IBP reduction. This means that the color- 
kinematic-dual representation of the four-loop form factor contains a truly 
free parameter! 

This detailed cancellation involves an order of fifty different master in¬ 
tegrals with up to ten propagators and occurs in both the planar and non- 
planar sectors. Most of the cancellations of the coefficients of the master 
integrals involve integrals out of all of the eighteen integral topologies which 
involve the free parameter and have nontrivial color factors. As an intricate 
example, consider the coefficient of the master integral 


INT[2, {-1,0,1,1,1,1,1, 0, 0, 0,1,1, 0, 0,1,0,0,0}] , 


(15) 


which occurs in topologies 12, 14, 15, 17 and 19. The coefficient in topology 
12 reads: 


-242 -h 2861e - 4906e2 - 37519e3 133706e^ - 185206^ - 378200e® -h 366000e^ 

21e4(l -K4e)(-1 -k5e)(l -F5e) 


(16) 

which is exactly —2 times the coefficient of this master integral in topology 
19. Taking into account the symmetry factors of topologies 12 and 19, 
which are 2 and 1, respectively, this leads to a cancellation. Similarly, the 
coefficients of this master integral in topology 14 cancels against that in 
topology 17, while the color factor of topology 15 is zero. The other master 
integrals typically involve much more intricate cancellations. 

It would be very interesting to reach a better understanding of the un¬ 
covered cancellations. On a practical level, this pattern completes the de¬ 
termination of the form factor integrand [1^. To our knowledge, this is the 
first example of a color-kinematic-dual integrand representation that con¬ 
tains a free parameter. It would be fascinating to see if there is any deeper 
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physical meaning of the free parameter, for example, as due to some gauge 
symmetries. This free parameter can be set to any value at the start of the 
computation. A natural value is —1, which simplifies the naive appearance 
of the integrals. This cancellation provides strong evidence for consistency 
of both the integrand result in and the IBP reduction performed in the 
present paper. 


4. Counting master integrals using algebraic techniques 


In 21[, Lee and Pomeransky developed a method to count the number 
of master integrals using a technique adapted from algebraic geometry. This 
is partially based on previous work by Baikov [^. This method only relies 
on the topology of the integral and, therefore, applies for arbitrary values 
of t and s. In particular, it is independent of obtaining an explicit IBP 
reduction and, hence, of the results of the previous section. This idea was 
also implemented in a Mathematica package called Mint 2^- In this section, 
we develop and apply this method to count the number of master integrals 
of four-loop form factors. 

Similarly to the start of IBP reduction, one first needs to determine the 


unique physical sectors for a given integral topology. As described in 2 ll |. 




it is convenient to perform this step using LiteRed [GOl. l6ll|. In particular, 
LiteRed determines the sector symmetries, which can be taken into account 
automatically in Mint. In the next step, the counting of master integrals is 
done independently for each physical sector. 

Let us consider a physical sector of L-loop topology and m propagators. 
From the associated polynomials U and F in equation ([9]), one can dehne a 
polynomial 

G{d) = U{a)+F{a), (17) 


new 


which has polynomial degree L -|- 1. Using the duality of homology and 
cohomology (see 211] for more details), the number of master integrals is 
related to the number of independent integral contours. The latter depends 
on the analytic structure of the polynomial G and is characterized by the 
so-called proper critical points. These are the stationary points at which G 


^Equivalently, one can consider a different polynomial using Baikov’s representation 
as discussed in [2l|, |^. For all cases being checked, the two polynomials always give 
the same result, while the use of the polynomial G is usually much more efficient in the 
computation. 
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does not vanish and are, thus, defined as the solutions of 


—— = 0 (z = 1 ,..., m) and G ^ 0. (18) 

UCXi 

Proper critical points may or may not exist for a given G polynomial. If 
there is no critical point, this implies that the corner integral of the sector 
is reducible, so that this sector does not contribute a master integral. On 
the other hand, if critical points do exist, the corner integral is irreducible 
with respect to IBP. This condition of irreducibility was first shown in [^ . 

A proper critical point may have a multiplicity which is measured by a 
topological invariant, the so-called Milnor number at the critical point. The 
number of master integrals is the sum over the Milnor numbers Mj of all 
proper critical points, 

^ of master integrals = Mi , (19) 

iGproper critical points 


provided that the proper critical points are isolated. We will explain the 
meaning of isolated critical point shortly. In the case of non-isolated critical 
points, more work is required 2]|, see also an explicit example in section 
KT\ below. 

Thus, the problem of counting master integrals is reduced to the one 
of computing the critical points of a system of nonlinear polynomial equa¬ 
tions given by equation (|18l) . which is usually highly nontrivial to solve. 
Fortunately, one can apply powerful algebraic approaches to make this job 
simple, based on the so-called Grobner basis technique. Readers who are 
not familiar with the Grobner basis and related concepts may consult 641 
for a pedagogical introduction. Here, only the basic procedure is outlined. 

The solution space of a system of polynomial equations is called affine 
variety associated with an ideal determined by the equations. For the system 
of equations (fT8]l . the ideal can be defined as 


/ = 


dG 

dai 


dG 
’ dam 


, aQG — 1 


( 20 ) 


where, in the last term, one introduces an additional parameter ao which 
forces the polynomial G to be non-vanishing at the critical point. The 
solutions of the affine variety can be obtained by computing the Grobner 
basis of equation ([20]) . 


gb(^) = {gi,g2,---,gk} ■ 


( 21 ) 
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Once the Grobner basis is obtained, it becomes relatively trivial to find the 
solutions. In our problem, we only concern ourselves with the number of 
solutions. This, in practice, can be conveniently obtained by counting the 
number of irreducible monomials in the Grobner basis]^ 

The above-described procedure is implemented in Mint 2ll| and applies 
straightforwardly to many simple cases; for example, applying it to the 
three-loop form factor, we obtain the result summarized in Table [2l The 
numbers of master integrals are consistent with the reduction of LiteRed. 


Table 2: Master integral counting for the six three-loop form factor topolo¬ 
gies given in Figure 5 of 13|. 


topology 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

# Mis via Mint 

9 

10 

10 

10 

14 

10 


For the four-loop form factor, however, one encounters two problems 
when trying to apply the Mint package. The first problem is that the set of 
critical points can form an affine variety of dimension > 1. This corresponds 
to the non-isolated critical-point case. Mint cannot deal with such cases 
automatically in its present version, but denotes them as Indeterminate. 
The number of such cases for the integral topologies at hand are summarized 
in Tabled Although the absolute number of these cases is low, they tend to 
occur in the more complicated sectors. They require some further work as 
described in [^. We employed a similar procedure and present an explicit 
example in subsection 14.11 


Table 3: Indeterminate sectors from Mint for the 34 topologies given in 
Figures 8 and 9 of 13]. 


topology 

(14) 

(15) 

(16) 

(17) 

(18) 

(20) 

(22) 

(23) 

(24) 

^ Indeterminate 

1 

1 

1 

1 

1 

2 

1 

1 

1 

topology 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(32) 

(33) 

(34) 

^ Indeterminate 

1 

3 

1 

1 

3 

1 

1 

1 

1 


Another problem is that, in quite a lot of cases. Mint gets stuck at the 
step of Constructing Grobner basis .... This is due to the complexity of four- 


®More precisely, one needs to compute the reduced Grobner basis, which is unique for 
a given monomial ordering, see e.g. [6J|. The number of irreducible monomials also takes 
the multiplicity of critical points (i.e. the Milnor number) into account. 
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loop integrals for which Mathematica cannot finish the computation of the 
corresponding Grobner basis. There are many known other packages which 
can compute a Grobner basis more efficiently. We applied Macaulay2 65l |. 
which can be conveniently used in Mathematica through mathematicaM2 .m 
(^ . and we employed Singular 671 in several most complicated cases. 

Based on Mint and with some extra effort to solve the above two prob¬ 
lems, we obtained the numbers of master integrals for the 34 different topolo¬ 
gies as summarized in Tabled 


Table 4: Master integral statistics: Mint results for the 34 topologies given 
in Figures 8 and 9 of jl3l |. 


topology 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

# Mis 

28 

35 

38 

34 

39 

48 

71 

52 

55 

69 

52 

95 

topology 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

# Mis 

76 

87 

95 

93 

106 

84 

105 

89 

45 

66 

41 

75 

topology 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 



# Mis 

55 

78 

92 

69 

93 

84 

35 

33 

31 

39 




Given the number of master integrals, Mint can also generate a set of 
master integrals. In most cases, a sector contains only one master integral, 
in which case the corner integral is chosen as the master integral. If there 
are more than one master integrals in a sector. Mint suggests integrals with 
higher powers of propagators (while Reduze typically chooses integrals with 
numerators). However, the choice of a double-propagator integral is rather 
heuristic; therefore, to find out whether it is a genuine master integral re¬ 
quires further checks. In our case, we could check them by using Reduze, as 
will be discussed shortly. 

To obtain the master integrals for the four-loop form factor, we still need 
to combine the master integrals of Table 0] together. Obviously, same master 
integrals may appear in different topologies of Table IH This may be detected 
by comparing their graph topology: if two graphs are isomorphic to each 
other, they are the same integral. After this step, we find that there are 
only 280 non-isomorphic master integrals: 244 with only simple propagators 
plus 36 containing a double propagator]^ 


®Note that we combined all 34 topologies including six zero-color-factor topologies: (8), 
(11), (15), (16), (18) and (20). If we exclude these six topologies, the number is reduced 
to 267, which is enough for the A/” = 4 SYM form factor. But for QCD with fundamental 
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Here, one should be cautious that there could be relations between inte¬ 
grals which are apparently not due to graph isomorphism, see e.g. 20, 5^ for 
some explicit examples. Such relations can be found by considering Feynman 
parameter representations of integrals and checking if the two representa¬ 
tions are related to each other by a permutation of Feynman parameters. A 
fast algorithm for such a check was proposed in 53|. Alternatively, one can 
also consider the so-called graph matroid technique to detect such hidden 
relations between different integrals, as implemented in Reduze 20(]. For us, 
there is an easier way to detect such relations, namely applying Reduze to 
reduce the basis integrals obtained by Mint. We checked that all 244 corner 
integrals correspond to basis integrals as determined by Reduze. Thus, they 
are indeed master integrals. 

As mentioned above, the choice of master integrals with double propaga¬ 
tors by Mint requires further checks. We applied Reduze to the 32 double¬ 
propagators integrals which do not appear in the twelve-propagator sector of 
topology 26 and found that four of the suggested master integrals are actu¬ 
ally reducible^ This does not change the number of Mint master integrals, 
but simply means that a different integral from this sector must be chosen 
as a master integral. By changing the position of the double propagator, we 
indeed found irreducible integrals. The twelve-propagator sector of integral 
topology 26 is analyzed further below. Thus, we obtained the hnal list of 
280 genuine master integrals. 

We would like to emphasize that the counting of master integrals based 
on Mint only relies on the topologies of given integrals and applies to ar¬ 
bitrary numerators; therefore, it applies to any theory, including QCD. We 
remind that the counting of Reduze given in Table [T] concerns particular 
master integrals appearing in the reduction of the A/" = 4 SYM form factor 
result. 


An interesting mismatch between Reduze and Mint 

Some caution is needed when assessing the output of Mint: this package 
is known to under-report the number of master integrals in very rare cases. 

To check this, we collected all master integrals appearing in the IBP 
reduction by Reduze, with t = 12 and s = 2 and with t = 13 and s = 1. It 
turns out that Reduze yields more master integrals. Many of them consist 


quarks, one likely needs all 280 master integrals. 

'^^Note that Reduze tends to choose integrals with numerators rather than donble prop¬ 
agators. So, irreducible means that, in the reduction resnlt of Reduze, there is an integral 
with the same propagators containing a numerator. 
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of twelve propagators and one to several numerators. They are at the edge 
of the reduction setting, and it is very likely that they may be reduced by 
including more IBP identities. These will not be discussed further here, 
details are available on request. What is unexpected is that there are nine 
corner integrals with only ten or less propagators, for example 

INT[2, {0,0,1,1,1,1,1,1,0,0, 0,1, 0,0,1,0,0,0}], (22) 

with only eight propagators. Mint shows that they contain no proper critical 
points, while Reduze and FIRE5 take them as master integrals. Due to the 
low number of propagators, most related IBP identities should already be 
taken into account by Reduze. Therefore, it is likely that these integrals 
are truly master integrals^ This would mean that the counting method of 
Mint requires some modihcation, although it is correct in the vast majority 
of cases. We cannot rule out, however, that there are integral relations 
which render this set of nine corner integrals reducible. It would be very 
interesting to understand this issue and see which is right. 

In addition, there is one more master integral in the Reduze basis, which 
appears to be an artifact of the way we solved the permutation symmetry 
consistency condition. 


4-1- Two worked-out examples 

Here, we give further details by considering two examples which cannot 
be computed directly by Mint. 

Let us first consider an example corresponding to an indeterminate case: 
a subsector of topology 26 with eleven propagators. 


INT[26, {1,1,1,0,1,1,1,1,1,1,1,1,0,0,0, 0, 0, 0}] . 


(23) 


From G = U + F, which are polynomials of {i = 1,..., 11) and of ho¬ 
mogeneity hve, one can construct the ideal I defined in equation (|20p . The 
corresponding Grobner basis cannot be computed by Mathematica, but can 
be calculated using Macaulay2 [6^. The quotient ^ace of the resulting 


ideal T = GB(/) turns out to be inhnite dimensionalo One then computes 
the primary decomposition of X (e.g. using Macaulay2), which gives 


X — X\ n To 


(24) 


^^Some other examples of two-loop integrals showing mismatch between Mint and IBP 
reduction were also found by V. Smirnov |68l| . 

^^The dimensionality of the quotient space is given by the number of irreducible mono¬ 
mials, which can be directly counted from the reduced Grobner basis as implemented in 
Mint [H. 


21 





where 


Xi =(186624ao — 3125, ai, 5a2 + 12, Sas — 12, 5a4 + 12,5a5 — 12, ag. Say + 12, 

5a8 + 12, Sag — 12,5aio + 12, 5aii + 12), 

X 2 =(729ao — 16, 2ai + 2a8 + 3, a 2 + 2aii + 6, a 3 — 3, 2a4 — 2a8 + 3, ag — 3, 

2ag “t" 2aii -t- 3, -t- 2a8 - 1 - 6,2aiia8 3a8 -{- 3a4i - 1 - 9,2a9 — 3,2aio — 2a4i -t- 3) . 


The quotient space of the hrst ideal is one-dimensional, which is obvious, 
since there is only one solution of the corresponding polynomial system. On 
the other hand, the quotient space of the second ideal is infinite dimensional. 
This can be treated as the example discussed in 21[. Rather than consider¬ 
ing the full ideal X2, one can simply consider the variety determined by the 
equation G(a8,aii) = 0, where 


(j(a8, an) — 2aiia8 -b 3a8 -b 3aii -b 9 


(25) 


is the only quadratic element in X2. Then, one computes the dimensionality 
of the quotient space of the ideal 

I = (^dG/das ,dG/dan , ohqG — 1'^ , (26) 

which turns out to be one. Altogether, one finds that the number of bases 
is two. 

As another illustrative example, we consider the sector of topology 26 
with full twelve propagators. In this case, the computation of the Grobner 
basis turns out to be hard even for Macaulay2. In this case, we apply 
Singular, which solves the problem in less than a minute using the method 
slimgb. The quotient space of the resulting Grobner basis I = GB(/) is 
seven-dimensional. One still needs to consider the symmetry which may 
reduce the number of independent master integrals. In order to make the 
symmetry obvious, one can compute the primary decomposition of I, which 
gives 

T = Xi n X 2 n X 3 n X 4 n X 5 n Xg, (27) 
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where 


Xi —(729o;o — 16, 4 q:i + 3, 0.2 — 3 ,03 — 3, 0 : 4 , 4 a 5 + 9, og — 3, otj — 3, 2q:8 + 9, Acng + 3, 
2q;io — 3, ctii + 6 , 2ai2 + 9), 

X 2 =(729ao — 16, cx\ — 3, 2 q ;2 + 9 ,03 — 3, 0 : 4 , + 6 , erg — 3, 4 cr 7 + 3, erg — 3, 2 er 9 + 9, 

2aio ~ 3,4erii + 9, 4 ai 2 + 3) , 

X 3 =(729ero — 16, 2eri + 9, er 2 — 3, 2 er 3 + 9, 2 er 4 — 3, erg + 6 ,4erg + 3, ery — 3,4er8 + 3, 
erg — 3, aio,4aii + 9, ai 2 — 3), 

^4 =^729erQ — 16, er^ — 3, 4 er 2 3, 4 er 3 3, 2 er 4 — 3, 4 er 5 4 “ 9,2erg 9, 2 er 7 9, erg — 3, 
erg — 3, aio, erii + 6 , eri 2 — 3), 

X 5 =(256ero + 1, eri + 2 , erg + 2 , erg + 2 , er 4 — 2 , erg + 6 , erg + 2 , ery + 2 , erg + 2 , ag + 2 , 
oio — 2 , erii + 6 , erig + 2 ) , 

Ig =(65536ero — 3125,5ai + 5 ai 2 — 8 , Sag + Soig — 8,03 — aig, 20 a 4 + Saig — 16, 
5ag + 16, ag — aig, Say + 5 ai 2 — 8 , Sag + Sayg — 8 , ag — aig, 

20 aio + 5ai2 — 16, 5aii + 16, 25a^2 “ 4:0ai2 + 128) . 

The quotient space of the last ideal is two-dimensional, while the others are 
all one-dimensional. So, naive counting would give seven master integrals. 

Topology 26, however, has permutation symmetry, see equation (HU). It is 
not hard to see that the hrst four ideals are all related to each other by this 
symmetry; therefore, they contribute only one independent master integral. 

In total, one Ends that there are four independent master integrals in this 
sector. For this integral topology, Reduze picks four twelve-propagator inte¬ 
grals as master integrals: the corner integral, two linear-numerator integrals 
and one with a quadratic numerator. Hence, the number of master integrals 
matches between Mint and Reduze. 


5. Towards numerical integration 


Given a set of master integrals, the next step is to perform the inte¬ 
gration. A signiheant part of the master integrals obtained in the previous 
sections are very complicated and seem to be on or slightly over the very 
edge of what current integration methods can achieve. There are various 
methods to perform numerical integration in the literature, see [l^ for an 
overview. 

It should be said that, in principle, Mellin-Barnes (MB) representations, 
which can be used very efficiently in solving planar topology integrals @,0 
Qizi, may be the optimal weapon of choice for four-loop integrals, if the 
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difficulties in obtaining efficient MB representations for non-planar integrals 
without hidden singularities can be resolved. See [l^, page 124, for a brief 
description of the problem. Progress in obtaining up to and including some 
three-loop non-planar representations in an automatic way was reported in 
[t^ ]. This may be used, at least in principle, to obtain MB representations 
for those integral topologies containing either a bubble or a triangle with 
two on-shell legs. 

We mostly explored numerical integration using sector decomposition 
In general, this is a method to turn a highly divergent integral into a 
sum over finite integrals, by resolving the singularities of the Feynman pa¬ 
rameter integral. These hnite integrals can then be integrated, for instance, 
using the CUBA library . We are aware of three different implementations 
of this idea in public codes: sector_decomposition [^, SecDec 77, 79| 

and FIESTA 81, 82|. Only the latter program in its latest versions (3.x) 
was able to handle some of the most difficult integrals of the problem at 
hand. 


5.1. Relevant implementation details 


There are various strategies for resolving the singularities of a given 
Feynman integral. FIESTA offers a choice of these. Most relevant here are 
strategies S (80(] and X [7^. Strategy KU (83|, which, in principle, is highly 
efficient in terms of sectors, did not resolve any of the sectors of the most 
complicated integrals (the corner integral of topology 25) in either of its 
three variants and was not pursued further. Strategy X, if it terminates, 
tends to yield better results in terms of the number of sectors than strategy 
5 . 


FIESTA splits the Feynman parameter integral in equation ([9]) into dif¬ 
ferent primary sectors. In each primary sector, one particular Feynman 
parameter is smaller than the others. If the integrand is invariant under a 
permutation of the a parameters, then these primary sectors may be related. 
A particular permutation symmetry of the a parameters descends from the 
parent graph. This symmetry may then be used to reduce the computa¬ 
tional workload. Let us recall the symmetry property of the scalar integral 
of topology 26 given in equation (|lip . From its orbits of the permutation 
group in equation (fT^ . 


{{1,2, 7,8}, {3,6,9,12}, {4,10}, {5,11}}. 


(28) 


we can see that, instead of computing primary sectors 1, 2, 7 and 8 sep¬ 
arately, one better just computes one of these, as the others will give the 
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same numerical answer. In practice, it can matter which coefficient is picked, 
since the structure of the integrand for the a parameters is different; this 
can result in different numbers of sectors. 

The integration time varies; for integrals with twelve propagators, the 
preparation time alone is typically measured in days to weeks. By contrast, 
all integrals with less than eleven propagators are highly tractable using 
FIESTA, typically requiring less than a day. 

Quadratic-numerator integral in topology 26 

The quadratic-numerator integral in topology 26 is naively expected to 
be one of the hardest integrals in the full set. It can be integrated using 
FIESTA. Integral topology 26 has two permutation symmetries, see equation 
(uni). To make use of this, one needs to pick a symmetric quadratic nu¬ 
merator. That is, one needs a numerator which is mapped to itself under 
both symmetries and which is related to the quadratic-numerator master 
integral in the output of Reduze by IBP reduction. Moreover, to facilitate 
the computation of the double derivative in the a parametrization, it is ad¬ 
vantageous to pick a numerator that is linear in the loop momenta. It is 
easy to check that 

{h-{pi-P2)f (29) 

transforms into itself under the symmetries in equation (jlOp . We checked 
that this is related to the quadratic master integral from the explicit IBP 
reduction. With the symmetries in place, one only needs to compute four out 
of twelve primary-sector coefficients in FIESTA. We picked the coefficients 
{4,0,4, 2, 2,0,0, 0, 0,0,0,0} for the order of propagators in equation (j7|l, 
basically choosing the lowest from the orbits in equation (j28p . 

The resulting integral can be integrated with some numerical effort, giv¬ 
ing the result 

INT26[(/3 • {pi - P2))^] =(-0.032986 ± 2 • lO”"^) e"® -b (0.0694456 ± 1.01 • 10"®) e"'^ 

-b (1.3506 ± 0.0002) e“® -b (-2.68804 ± 0.00317) e"® 

-b (-6.23707 ± 0.04013) -b (12.6763 ± 2.0782) e"^ 

-b (1234.49 ± 32.97) e"VO(e"^). (30) 

Preparing for the numerical integration takes about five days. The inte¬ 
gration time itself strongly depends on the integration parameters in the 
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Vegas integration algorithm of the CUBA library^ We checked by running 
with different numbers of evaluation points that the first five coefficients are 
numerically stable under variations of the relevant parameters, especially 
maxeval. In other words, it is likely that their errors are in the reported 
ranges. 


6. Summary and discussion 


In this paper, we obtained an explicit integral reduction for the four- 
loop Sudakov form factor in A/" = 4 SYM based on the integrand result in 
|l3| | and by exploiting the Reduze package [^. This reduces the challenging 
four-loop problem, in particular the computation of the unknown non-planar 
cusp anomalous dimension, down to the evaluation of an explicit set of given 
master integrals. 

We also obtained master integrals by using the algebraic method intro¬ 
duced in 21[ . This provides a first non-trivial implementation of the method 
to a complete observable. Limits of the Mint package and their resolutions 
are explained. The resulting integrals are, with some qualifications, the 
master integrals of quite generic form factors. For instance, they are ex¬ 
pected to cover all integrals which appear in the QCD result. By comparing 
with the results of IBP reduction, we found interesting mismatches, namely, 
existing masters integrals from IBP reduction are taken to be reducible by 
the Mint package. It would be important to understand the reason of this 
discrepancy. 

Based on the reduction result, we furthermore observed surprising can¬ 
cellations in the form factor. First, the form factor integrand obtained in 
[l^ contains a free parameter. While it was unclear whether this parame¬ 
ter could be uniquely fixed by further physical constraints, we showed that 
the terms depending on this parameter all magically cancel upon using IBP 
identities. Therefore, this parameter is truly free and gives the first example 
of a color-kinematic-dual integrand with such a free parameter. Further¬ 
more, two integrals, corresponding to topologies 18 and 20, turn out to be 
exactly zero, which is not obvious at all at the integrand level. It would be 
very interesting to have a physical understanding of these intriguing cancel¬ 
lations. 

What remains to be done is to compute the master integrals. However, 
for many of them, this is currently just out of reach for publicly available 


brief exploration of the other algorithms in the CUBA library quickly shows that 
Vegas is the optimal choice here, as for both time and accuracy. 
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(implementations of) integration methods. There is an important analogue 
to the three-loop case, which we would like to draw attention to. In that 
case, a reduction to master integrals was obtained in 0]. Their integration 
was first performed almost three years later in [^, using partly numerical 
methods. An analytic result was then reported in 8^, with an important 


cross-check in 8^. Currently, the three-loop form factor integrals often ap¬ 


pear as an example application of new integration and organization methods, 
e.g. in 86, 87|. 


A new feature in the four-loop case is the appearance of several master 
integrals beyond corner integrals. These can be taken to be either doubled- 
up propagators or single numerators. The latter can, in principle, be traded 
for higher-dimensional scalar integrals using the Feynman parametrization 
in equation ([9j) . The problem is that these scalar integrals feature four dots: 
there is a sum over integrals with four doubled-up propagators. Without 
an IBP reduction of these integrals down to the basis, this is an unpractical 
result. An explicit IBP reduction with up to four dots would enable the use 
of several methods. For instance, one could use dimensional recurrences 
to either compute the integrals numerically or to express the form factor 
in terms of a quasi-finite basis 87| . An explicit quasi-hnite integral reduction 


for the three-loop form factor has been announced 90|. The needed IBP 


reductions at four loops seem very hard to obtain within the current version 
of Reduze and simply out of reach for anything else on the public market. See 


91[ for a proposal of an improved algorithm which is particularly suitable 


for single-scale reduction problems. 

There is a practical issue when calculating the form factor which concerns 
the choice of master integrals. There are several choices of master integrals 
one could make as for both linear numerators and doubled propagators. 
Which one to choose should then be guided by integrability, i.e. by which 
of these choices they can actually be integrated. Moreover, as for numerical 
evaluation, there are issues of precision and speed to be considered. One 
would expect doubled-up propagators to give simpler Feynman parameter 
integrals, while linear numerators tend to give less divergent coefficients. A 
canonical example of the two choices are the master integrals of the double¬ 
box-integral topology. One can either choose a doubled-propagator basis 
element, as first done in 9^. Alternatively, one can use a linear-numerator 


integral, as in 93|]. The latter has a less divergent prefactor as a function 


of e. Which choice is the more practical must be decided on a case-by-case 
basis though; we observed cases where no difference in expansion order is 
present. 

There are various other avenues to be explored which may lead to a 


27 












full computation of the cusp anomalous dimension through the form factor 
route. One would be to explore MB representations of four-loop non-planar 
integrals. An obstacle here is that naive implementations of this idea quickly 
lead to high-order MB representations containing an order of one hundred 
MB integrations. Any algorithm which could curb this number would offer a 
way forward. While numerical evaluation may have the precision issue, ana¬ 
lytic methods would be more desirable. Linearly reducible integrals may be 
evaluated analytically by Hypint 941]. This works particularly efficiently for 


integrals of less singular behavior 87|. An alternative route for computing 


the cusp anomalous dimension is to extend integration techniques based on 
differential equations 95|] to form factors by making one extra leg off-shell 
and then taking the massless limit. Controlling this limit after integration 
might also be an obstacle. Apart from that, IBP reduction, a necessary 
step in constructing the differential equation, would be more challenging for 
deformed integrals. See [86| for an application to the three-loop case. 
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Appendix A. Three-loop form factor check 

The three-loop form factor was first computed in [^; we started from 
the expression obtained in [I^, which can be shown to be equivalent. The 
problem consists of six integrals with at most a single numerator. The 
integrals contain nine propagators, and three additional numerators have to 
be specihed to obtain a dehnite basis. Five of these integrals have a single 
permutation symmetry. 

IBP reduction using Reduze on a part of our local computing cluster 
running in parallel mode took about 1.5 days. The reduction yielded a set 
of master integrals, the most complicated of which indeed correspond to 
those obtained in [^. These integrals were then computed using FIESTA, 
running for a few days with their standard precision setting. The result for 
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the form factor is 


pi^) _ 

2,num 


0.166606 


0.00012 0.00054 1.10182 0.93788 ^/l\ 

(A.l) 


which is to be compared to the known result, e.g. from 


i, 


7(3) 

2, exact 


1 , 11C(3) , 2477r4 

12e3 25920 


+ 0 



(A.2) 


The difference between the exact answer and the numerical approximation 
is 


(3) 0.000060 0.00012 0.00054 0.00006 0.00964 „/1\ 

^^2 =— - — + — -+ + 

(A.3) 

Note that the final coefficient is only off by about one percent. With more 
numerical effort, this could be brought down more. For the purposes of 
this article, it is important that the calculational pipeline works, at least in 
principle. 


Appendix B. Alternative choice of numerators 

The following choice of numerator factors for master integrals optimizes 
the power series expansions of the coefficients to some extent. This choice 
is guided by the following principles: 

• choose the numerators linear in the loop momenta to minimize the 
complexity of the Feynman parameter integrals; 

• when expressed in the master integral basis preferred by Reduze, the 
coefficient of the linear-numerator master integral comes with a nega¬ 
tive power of e if possible; 

• when expressed in the master integral basis preferred by Reduze, the 
coefficient of the scalar version of the linear-numerator master integral 
comes with a negative power of e; 

• preserve permutation symmetry if possible. 
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This leads to the following list, where the coordinates are in the basis of 
the integral topologies given in appendix Appendix C[ 

INT[30, {1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0, 0, 0},p2 • k] , 

INT[3, {1,0,0,1,0,1,1,0,1,1,0,1,0,1,1,1,1, l},pi -Ga-Zs)], 

INT[28, {1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0, 0, 0},p2 • ^ 5 ], 

INT[26, {1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0, 0, 0},pi • h] , 

INT[25, {1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0, 0, 0}, {p2 - Pi) ■ h] , 

INT[23, {1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0, 0, 0},p2 • ^ 5 ], 

INT[20, (1,0,1,1,1,1,1,1,1,1,1,1,0,0,1,0, 0, 0},p2 • ik + k- /s)], 
INT[17, (1,0,1,1,1,1,1,1,1,1,1,1,0,0,0,1, 0, 0},p2 • ih - k)] , 

INT[11, (1,0,1,1,1,1,1,1,0,1,1,1,0,0,1,0, 0,1},Pi -/s], 

INT[5, (1,0,1,1,1,1,1,1,1,1,1,1,0,0,1,0,0,0}, (p2 - Pl) • h] , 

INT[4, (1,0,1,1,1,1,1,1,1,1,0,1,1,1,0,0,0,0},pi -y, 

INT[3,{1,1,0,1,1,1,1,0,1,1,0,1,0,1,1,0,1,0},P2 • {k - k)] , 
INT[2,{0,0,1,1,1,1,1,1,1,1,1,1,0,0,1, l,0,0},pi ■ k] . 


This choice leads to a signihcant reduction in the IBP computational work¬ 
load, at the expense of introducing non-symmetric numerators. 


Appendix C. Basis 

This appendix contains the main basis used in the reduction. The num¬ 
bering of the equations corresponds to the topologies in [l^. In each case, 
the twelve entries given in the first line parametrize the twelve propagators 
of the respective integral and the six entries in the second line the chosen 
numerators. Note that the numerators do not respect the symmetries of the 
diagram topology corresponding to the first twelve propagators, with the 
exception of topology 1. We have dehned q = Pi + P2- 


{k, ki k, h, k — k-, k — h-, Pi — k-, k — k, —k + q-, —k + q, —h + q-, —k + q, 
h — k, k — k, k — k, k — P2, k — P2, k — P2} , (C.l) 

1 ^ 6 ) k, k, k, k—k, k—k, Pi—k, k—k, —k+q, —k+q, —k+q, —k+k—k+q 
k — k, k — k, k — k, k — Pi, k — P2, k — P2} , (C.2) 

1 ^ 6 ) k, k, k, k—k, k—k, k—k, Pi—k, —k+q, —k+q, —k+q, —k+k—k+q 
k — Pl, k — P 2 , k — P 2 , k — P 2 , k — k, k — k} , (C.3) 
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{^6; h, hi Pi—hi h—hi h—hi h—hi —h+h+P2i —h+Qi —h+Qi —h+Qi 
h — hi h — hi U — hi h — P 2 i h — P 2 i h — P 2 } , (C.4) 

{^6; hi hi hi h—hi h—hi Pi—hi h—hi —h+Qi —h+h+P2i —h+Qi —h+h—h+Qi 
h — P 2 i h — hi h — hi h— Ph h— P 21 h — P 2 } , (C.5) 

{^6; hi hi hi Pi—hi h—hi h—hi P2—hi —h+Qi —h+Qi —h—h+Qi —h—h+Qi 
h — hi h — hi h — hi h — hi h — P 21 h — P 2 } , (C.6) 

{^6; hi hi hi Pi—hi h—hi h—hi P2—hi —h+Qi —h+h+hi —h+Qi —h—h+Qi 
h — hi h — hi h — hi h — P 21 h — P 21 h — P 2 } , (C.7) 

{^6; hi hi hi P2—hi h—hi h—hi Pi—hi -h+Qi —h+Qi —h—h+Qi h—h+h—P2i 
h — hi h — hi h — hi h — hi h — P 21 h — P 2 } , (C.8) 

{^6, hi hi hi h—hi Pi—hi h—hi —h+Qi —h+Qi —h—h+Qi h+h—Qi —h—h+Qi 
h — hi h — P 2 i h — P 2 i h — hi h — hi h — P 2 } , (C.9) 

{^6, hi hi hi h—hi Pi—hi h—hi h—h+hi —h+Qi —h+Qi —h—h+Qi h—h+h+h—Qi 
h — hi P2 — hi h — hi h — P2i h — P2i h — h} , (C.IO) 

{^6; hi hi hi h—hi Pi—hi h—hi —h+Qi —h+Qi h+h—Qi —h—h+Qi h—h+h+h—Qi 
h — hi h — P 2 i h — Ph h — P 2 i h — hi h — P 2 } 1 (C.ll) 

{^6; hi hi hi P2—hi Pi—hi h—hi h—hi —h+Qi —h+h+hi —h—h+Qi —h+h—h—h+Qi 
h - hi h - P2i h - P2i h - Ph h - P 21 h - Pi} , (C.12) 

{^6; hi hi hi h—hi P2—hi Pi—hi h—hi —h+Qi —h—h+Qi —h—h+Qi —h—h+Qi 
h - hi h - P2i h - P2i h - Ph h - P 21 h - Pi} , (C.13) 

{^6; hi hi hi P2—hi Pi—hi h—hi h—hi —h+h+hi —h+Qi —h—h+Qi —h—h+Qi 
h — hi h — hi h — hi h — P 21 h — P 21 h — P 2 } , (C.14) 

{^6; hi hi hi P2—hi Pi—hi h—hi h—hi —h+Qi —h—h+h+Ph —h—h+Qi —h—h+Qi 
h - hi h - P2i h - P2i h - Ph h - P 21 h - Pi} , (C.15) 

{^6; hi hi hi P2—hi Pi—hi h—hi h—hi —h+Qi h—h+h—hi —h—h+Qi —h—h+Qi 
h — hi h — P2i h — hi h — P 21 h — Ph h — Pi} , (C.16) 

{^6; hi hi hi P2—hi h—hi Pi—hi h—hi —h+Qi —h—h+Qi —h+h—h+Qi —h+h—h—h+Qi 
h - P 2 i h - hi h - P 2 i h - Ph h - P 2 , h - Pi} , (C.17) 

{^6; hi hi hi Pi—hi P2—hi h—hi h—hi h—h—hi —h+Qi —h+h—h+Qi —h+h—h+h+Ph 
h — P2i h — hi h — hi h — P 21 h — P 21 h — Pi} > (C.18) 
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{^6; h, h, hi P2—hi Pi—hi h—hi h—hi —h+h+hi —h+Qi —h+h+hi —h—h+Qi 
h — hi h — hi h — Ph h — P 2 i h — P 2 i h — P 2 } , (C.19) 

{^6; hi hi hi P2—hi h—hi Pi—hi h—hi —h+h+hi —h+Qi —h+h+hi —h—h+(li 
h — hi h — hi h — Ph h — P 21 h — P 21 h — P 2 } , (C.20) 

{^6; hi hi hi Pi—hi h—hi P2—hi h+hi —h—h+Ph —h+Qi —h+h+Ph —h—h+Ph 
h — P2i h — P2i h — P2i h — Ph h — hi h — h} 1 (C.21) 

{^6; hi hi hi h—hi h+hi P2—hi Pi—hi —h+h+Ph —h+Qi —h—h+Ph —h+h+h+hi 
h — hi h — P2i h — P2i h — Ph h — P 21 h — Pi] 1 (C.22) 

{^6; hi hi hi h—P2i Pi—hi h—hi h—hi —h+Qi h—h—hi —h+h+Ph —h+h+Ph 
h — P2i h — P2i h — Ph h — Ph h — P2i h — h] 1 (C.23) 

{^6, hi hi hi P2—hi h+hi h—hi pi—hi —h—h+Ph —h—h+Ph —h+Qi —h—h—h+Qi 
h — P2i h — hi h — P2i h — Ph h — P 21 h — Pi] , (C.24) 

{^6, hi hi hi Pi—hi h—hi P2—hi h—hi —h+h+hi —h+Qi —h+h—h+Ph —h+h—h—h+Qi 
h — hi h — hi h — hi h — Ph h — P 2 i h — P 2 ] , (C.25) 

{^6; hi hi hi h—hi h—hi Pi—hi P2—hi —h+h+hi —h+Qi h—h+h—P2i 

— h + h — h — h + Qi h — hi h — hi h — hi h — Ph h — P21 h — P2} 1 

(C.26) 

{^6; hi hi hi h+hi P2—hi Pi—hi h—hi —h—h+Ph —h+Qi —h+h+h+P2i —h+h+h+hi 
h — hi h — P2i h — P2i h — Ph h — P 21 h — Pi] 1 (C.27) 

{^6; hi hi hi h—hi Pi—hi P2—hi h—hi —h+Qi h—h+h—P2i —h—h+Qi —h+h—h+P2i 
h — P2i h — hi h — hi h — P 21 h — P 21 h — Pi] 1 (C.28) 

{^6; hi hi hi Pi—hi h—hi P2—hi h+hi —h—h+Ph —h+Qi —h+h—h+P2i —h—h—h+Qi 
h — P2i h — hi h — P2i h — Ph h — P 21 h — Pi] 1 (C.29) 

{^6; hi hi hi h—hi h+hi P2—hi Pi—hi —h+Qi —h—h+Ph —h+h+h+P2i —h+h—h+P2i 
h - P2i h - hi h - P2i h - Ph h - P 21 h - Pi] , (C.30) 

{^6; hi hi hi h—Ph P2—hi h+hi h—h—P2i —h—h+Ph —h+Qi h+h—P2i —h—h+Qi 
h - hi h -Ph h-P2i h -Ph h - P 21 h - Pi] , (C.31) 

{^6; hi hi hi h—Ph P2—hi h+hi —h—h+Ph —h+h+P2i —h+Qi h—h—hi —h—h+Qi 
h - hi h -Ph h-P2i h -Ph h - P 21 h - Pi] , (C.32) 

{^6; hi hi hi P2—hi h+hi h—Ph —h+Qi h+h—P2i —h—h+Qi h+h+h—P2i —h—h—h—h+Qi 
h - hi h -Ph h-P2i h -Ph h - P 21 h - Pi] , (C.33) 
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{^6) hi h, h, P2—h, h—pi, h+le, h+h—P2, —h—h+Pi^ —h+Q: —h—h+q, —h—h—h—h+Q: 
h — ^6) 4 — Pi, h — P2, h — Pi; h — P2, h — Pi} ■ (C.34) 

Appendix D. List of ancillary files 

• List of master integrals appearing in the form factor obtained by 
Reduze 

• List of master integrals as suggested by Mint 

• IBP reduction of all integrals appearing in integral topology 25 for the 
A/" = 4 SYM form factor 
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